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ABSTRACT 

Aims. We investigate the properties of the Galactic ISM by applying Faraday tomography to a radio polarization data set in the 
direction of the Galactic anti-centre. 

Methods. We address the problem of missing large-scale structure in our data, and show that this does not play an important role for 
the results we present. 

Results. The main peak of the Faraday depth spectra in our data set is not measurably resolved for about 8% of the lines of sight. An 
unresolved peak indicates a separation between the regions with Faraday rotation and synchrotron emission. However, cosmic rays 
pervade the ISM, and synchrotron emission would therefore also be produced where there is Faraday rotation. We suggest that the 
orientation of the magnetic field can separate the two effects. By modelling the thermal electron contribution to the Faraday depth, 
we map the strength of the magnetic field component along the line of sight. Polarized point sources in our data set have rotation 
measures that are comparable to the Faraday depths of the diffuse emission in our data. Our Faraday depth maps show narrow canals 
of low polarized intensity. We conclude that depolarization over the telescope beam produces at least some of these canals. Finally, 
we investigate the properties of one conspicuous region in this data set and argue that it is created by a decrease in line-of-sight 
depolarization compared to its surroundings. 

Key words. Magnetic fields - Polarization - ISM: Magnetic fields - Radio continuum: ISM - Techniques: polarimetric 



1. Introduction 

Faraday rotation, the rotation of the plane of linear polarization 
due to the birefringence of a magneto-ionic medium, provides a 
powerful tool for exploring the magnetic universe, complement- 
ing effects such as Zeeman splitting and the alignment of ellip- 
soidal dust grains perpendicular to magnetic field lines. From 
external galaxies we have learnt about the global structure of the 
magnetic field, about its alignment with the galactic spiral arms, 
and of the importance of the magnetic pressure in ISM dynam- 
ics, which is comparable to the thermal pressure in the cold and 
warm neutral phases of the ISM (see Beck 2007a, 2007b). 

In the Milky Way, rotation measures of pulsars, combined 
with dispersion measures, yield a rather complicated picture of 
the large-scale Galactic magnetic field, with evidence for sev- 
eral field reversals (see e.g. Han et al. 2006 ). Rotation measures 
of extragalactic sources from the International Galactic Plane 
Survey (IGPS) indicate a simpler structure of the large-scale 
field, without much evidence for many field reversals (Brown 
et al. l2007l i. 

Studies of Faraday rotation of the diffuse Galactic ISM offer 
a big advantage over other investigations based on pulsars and 
extragalactic sources, which is that the diffuse emission is visible 
in all directions. This has been spectacularly demonstrated by 
the 1.41 GHz survey carried out with the DRAO 26 m telescope 
by Wolleben et al. ( I2006I I. A milestone in this type of work is 
the polarization survey by Brouw & Spoelstra in the late 1970s 
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( 119761 and Spoelstra 1 1984| l. which has the largest sky coverage 
in rotation measures of the diffuse emission to date. Since the 
1990s interferometric studies at low frequencies have led to new 
insights in e.g. the properties of the ISM on small scales, and on 
turbulence in the Galactic ISM (Wieringa et al. ,1993; Haverkorn 
et al. 120041 |2006t . 

Recently, a method called Faraday tomography has been de- 
veloped with which one can study the relative distribution of 
synchrotron-emitting and Faraday -rotating regions along the line 
of sight (see Brentjens & De Bruyn 2005 ). Application of this 
technique for studying the Galactic ISM is a very new field of 
research, and not much literature is available on this subject. In 
Schnitzeler et al. (2007b) we discussed a low-frequency data set 
that covers » 50° in the direction of (/,/?) = (181°, 20°), that we 
obtained with the WSRT. We concluded that many lines of sight 
contain at least one Faraday screen, a region in which Faraday 
rotation and synchrotron emission are not mixed. This result 
appears counterintuitive, since one would expect synchrotron- 
emitting cosmic rays to be present everywhere there are mag- 
netic fields. In this paper we present the full analysis of this data 
set. 

Faraday tomography promises to become an interesting and 
exciting way to investigate the Galactic ISM that will become 
more and more available in the near future, with the advent of 
new surveys and radio telescopes that combine high spatial res- 
olution with high frequency resolution, like the GALFACTS sur- 
vey that is carried out with the Arecibo telescope (Taylor 2004*, 
the LOFAR aiTay (Rottgering et al. ,2006i ). and the MWA SKA- 
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precursor (Bowman et al. |2006l l. A better understanding of the 
Galactic foreground will also be useful for studies of the cosmic 
microwave background. 

This paper is the first in a series of articles in which we dis- 
cuss different regions in the second Galactic quadrant. Here we 
discuss some of the observational features that we encounter, 
and we extend some of the techniques that have been used to 
study the diffuse Galactic ISM with a relatively small number 
of frequency channels to the many-channel regime. First, we 
give a short overview of Faraday tomography in Sect. |2] and 
we illustrate its potential with a couple of examples. In Sect. 
[3] we discuss the observational characteristics of the GEMINI 
data set, and in Sect. |4] we present the GEMINI data. We show 
that Faraday modulation in the foreground can convert polarized 
emission on large angular scales to smaller angular scales that 
can be picked up by an interferometer that has no short baselines. 
We derive the magnetic field strengths implied by the observed 
Faraday depths of the diffuse emission in Sect. |5] and in Sect. |6] 
we compare the Faraday depths we derive for the GEMINI data 
set with the rotation measures that we find for polarized point 
sources. One conspicuous feature in the GEMINI data are the 
dark canals in the maps at constant Faraday depth, and we inves- 
tigate their origin in Sect.|7] In Sect. [8] we present an explanation 
for a large and bright polarized region in our data set. 



2. Notes on Faraday tomography 

2. 1 . Faraday tomography 

The wide coverage in space of our data (from 0.6 m^ < 
< 0.9 m^), combined with the relatively small channelwidth 
{6A^ 1.5x10"^ m^), allows us to do Faraday tomography, also 
known as Rotation Measure Synthesis (see e.g. Brentjens & De 
Bruvn l2005T l. In this section we introduce the concept of Faraday 
tomography, and we illustrate how it can be used to study the 
properties and the relative distribution of regions with Faraday 
rotation and synchrotron emission along the line of sight. 

The observed Stokes parameters Q(A) and U(A) are the 2 
orthogonal components of the polarization vector P{A) - QiA) + 
iU{A). P(A) is the vector sum of all polarization vectors that are 
emitted along the line of sight and that are Faraday rotated in the 
ISM between the point of emission and the observer: 



00 00 
P(A) = j P(x) e^i'^W-^'dx = j i'CR)e''^''' 



(1) 



where the first integral is over physical distance 'x', and the 
Faraday depth of point 'x' 



nix) [rad/m^] = 0.81 



bserver 



[cm-^] B IjiG] ■ dl [pc] (2) 



source at 'x' 



measures the total amount of Faraday rotation between the point 
of emission 'x' and the observer. In the second integral of Eqn.[T] 
we have replaced the integral over physical distance by an inte- 
gral over Faraday depth. If there is a one-to-one correspondence 
between 'x' and "R, this change in coordinate is trivial. However, 
a change in the direction of the magnetic field component along 
the line of sight B|| = B ■ dl assigns the same fl to different 
'x'. To go back from fl to 'x' thus requires assumptions on the 
magnetic field geometry. In a future article we will return to this 



issue. Faraday tomography is based on the inversion of the sec- 
ond integral in Eqn.[T] 



PCR) ^ K J W(A^) P(A^)s^^'''^'''dA^ 
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where 




W(A^) dA^ 



(3) 
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produces the correct normalization for PCR). W(A^) = 1 if the 
A^ has been observed, otherwise W(A^) = 0. This inversion can 
be interpreted either as a 'coherent addition' of the observed 
P(A) by using an assumed Faraday depth to derotate the po- 
larization vectors, or, similarly, as a Fourier transform from A^ 
space to K space. Note that the integral is also calculated over 
negative values of A^, which is physically impossible. If the 
measured P{A^) are in K (or mJy/beam), then also the PCR) 
will be in K (or mJy/beam), or, more precisely, in K/RMSF 
width (or mJy/beam/RMSF width). If we want to compare dif- 
ferent data sets, with different RMSFs, it is necessary to convert 
this quantity to K/rad/m^ instead of K/RMSF (equivalently for 
mJy/rad/m^). Here we only consider one data set, and we did 
not convert the measured PCR) to K/rad/m^. 

The complex polarization vector P(A) can be written in terms 
of the observed polarized intensity P(A) and polarization angle 
<D(^) as P(A) = P(A) e^'^'^'^'^K Similarly, PCR) = PCR) e^'**'^', 
where PCR) is the intensity of the polarized emission at Faraday 
depth K, and <t>CR) is the orientation of the electric field vector of 
the synchrotron radiation emitted at Faraday depth K, which is 
perpendicular to the local direction of the magnetic field at that 
Faraday depth. A Faraday depth spectrum (or "R spectrum) can 
be constructed by calculating PCR) and OCR) for many "R values. 

The Fourier-transform nature of Eqn.[3]implies that it shares 
some characteristics with other Fourier-transform based meth- 
ods like radio interferometry. Eqns. 61 to 63 from Brentjens & 
De Bruyn (120051 ) describe this behaviour quantitatively, and for 
convenience we reproduce these equations here. We assume here 
that the observations uniformly span the A^ range from A^^^ to 
AL.,^ - A- ■ -I- A/l^ with weight 1 , and with weight outside this 

max mm ^ ' ^ 

range. 

The finite extent in the AA^ coverage of the data introduces 
an instrumental response along the axis, known as the RMSF 
(Rotation Measure Spread Function), with which the PCR) spec- 
trum is convolved. The shape of the (normalised) RMSF is given 
by 



RMSF(??) = K I WiA^)e-^''^-^^dA^ 



(5) 



. The FWHM of the RMSF depends on the total range in A^ space 
covered by the observations, AA^, according to 



3 8 

FWHM [radians /m^] = 

AA^ 



(6) 



. We replaced the factor of 2 V3 from the original paper by 
Brentjens & De Bruyn by 3.8, because 3.8 is closer to the 
FWHM of the sine RMSF response than 2 V3 -3.5. 
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Fig. 1. P(A^) and <S>(A'^) distributions (top 2 panels), and PC??) and 
<i>('R.) spectra (bottom 2 panels) for a line of sight with a barely 
resolved main peak that is similar to a Faraday screen. In the 
text we describe how we derotated the polarization angles in the 
bottom panel of this figure. The x^^^ of the RM fit is 0.99 (de- 
termined by using the RM with the maximum f (RM)). The K 
spectra are sampled from -100 rad/m^ to -i-lOO rad/m^ in steps 
of I rad/m^. In the bottom 2 panels, soUd lines indicate the cal- 
culated PCR) and OCR) spectra, and the dashed lines indicate the 
fitted RMSF. 

In general the observations will not cover all wavelengths 
down to m^, which means that "R-extended structures will be 
missing from the PCR) spectra. This is similar to the missing 
large-scale structure problem interferometers suffer from when 
observing extended emission. If the source has a Gaussian PCR) 
distribution along the line of sight, Brentjens & De Bruyn cal- 
culate that in the reconstructed PCR) spectrum "R scales larger 
than 



maximum scale [radians /m ] ^ — — 

A ■ 

mill 



(7) 



are suppressed by more than 50% by the missing short wave- 
lengths. 

Analogous to how in an interferometer the size of the indi- 
vidual dishes sets the field of view, the channelwidth expressed 
in A^, 6A^, determines the K for which the sensitivity has dropped 
to 50% due to smearing of the polarization angles over individ- 
ual channels, "Rmax- 



7 1.9 

I'^maxl [radians/m ] = — 



(8) 



. Also in this equation we replaced the V3 from the paper by 
Brentjens & De Bruyn by 1.9 because this is closer to the cor- 
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Fig. 2. Identical to Fig. [T] but for a line of sight towards a po- 
larized extragalactic source that also contains diffuse Galactic 
emission. The;^'^^^ of the RM fit is 8.72. 



rect value for the HWHM of the sine RMSF response. Brentjens 
& De Bruyn discuss the (minor) effect of using (possibly non- 
contiguous) channels with constant width in frequency instead 
ofA^. 

The addition of polarization vectors originating at different 
distances along the line of sight in general leads to depolariza- 
tion, either because the vectors are emitted with different orien- 
tations and/or because different amounts of Faraday rotation in- 
duce misalignment of the polarization vectors. Faraday tomog- 
raphy can separate the emission coming from subregions with 
different H along the line of sight, thereby reducing the influ- 
ence of these sources of depolarization. 

Brentjens & De Bruyn (2005) note that one can also derotate 
the polarization angles of PCR) to A^ + 0. In that case the A^ 
in the complex exponentials of Eqns. [3] and |5] are replaced by 
{A^ - /Iq). Brentjens & De Bruyn showed that derotating to the 



weighted average of the observed A , 



J W{A^)A^dA^ 



J W(A^)dA^ 



(9) 



minimises the polarization angle variation over the main peak of 
the RMSF. However, now the polarization angle at the peak of 
the RMSF at Faraday depth !Rpeak is no longer perpendicular to 
the magnetic field component in the plane of the sky. To deter- 
mine the orientation of this magnetic field component, one needs 
to correct the polarization angles O by 
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OoCR) = ^(m-'Rp..^^ (10) 

. We derotated the polarization angles in the bottom panels of 
Figs.[T]and|2]in this way. 

2.2. Examples 

We illustrate for a number of geometries of synchrotron-emitting 
and Faraday -rotating regions the resulting P(A) and as well 
as the PCR) that is reconstructed from the observations. 

First, consider the case where there is only one infinitely thin 
source of emission along the line of sight, at Faraday depth : 
PCR) - Pi dCR - Hi), where 6{x) is the Dirac delta function. 
Then P{A) = Pi e^'^"*': the P{A) are identical for all observing 
wavelengths A, where we assumed that the synchrotron emission 
has a spectral index of 0. As Brentjens & De Bruyn (120051) have 
discussed, a non-zero spectral index introduces only a distortion 
in the RMSF away from the main peak. O depends linearly on 
A^, and since the rotation measure RM = d 0/5 A^, RM = H 
in this case. The Pi'R) spectrum that is reconstructed from the 
observed P{A) and ^{A) is identical to the original Pi'R) spec- 
trum, convolved with the RMSF. The situation described in this 
paragraph occurs in our data: in Fig. [1] we show the measured 
P{A^) and reconstructed Pi'R) spectra for a line of sight with an 
essentially unresolved Pi'R) spectrum. 

Next, consider the case where there are two peaks of height 
P\ and P2 at Faraday depths 'R 1 and 'R 2. In this case PiA) - 
p^g2iR i/i-_|.p^g2i« 2-1- jjjjg configuration produces a beat in PiA), 
because the 2 polarization vectors of length Pi and f 2 rotate at 
different speeds 'R 1 and 2 in the Stokes iQ,U) plane. Note that 
cD no longer depends linearly on A^, which means that RM will 
not be constant. The Pi'R) spectrum that is reconstructed from 
the observed PiA) will show the sum of 2 ^-functions convolved 
with the RMSF. The peaks in the reconstructed Pi'R) spectrum 
will not lie at exactly 'R i and "R 2, because the presence of the 
other peak will influence the shape of both peaks. Therefore, 
since <t>(:R) is derived for the Faraday depth of the peak in the 
reconstructed Pi'R) spectrum, the derived <i>i'R) will not be the 
intrinsic position angle of the electric field of the emitted syn- 
chrotron radiation. When a polarized extragalactic point source 
is producing one of the 2 peaks, and the diffuse Galactic emis- 
sion produces the second peak, the x^^^i of the linear fit of <1> 
versus A^ will be very high. In Fig.|2]we illustrate this case. The 
peak afR = -28 rad/m^ is produced by an extragalactic source, 
and the region around 'R - +5 rad/m^ is produced by Galactic 
emission in the foreground, as we established by comparing this 
line of sight to adjacent lines of sight, where the extragalactic 
source is not present. 

The previous two examples dealt with infinitely thin emis- 
sion regions. If the emission region has a finite depth, but is 
not mixed with Faraday-rotating ISM, the results we discussed 
in the previous sections still hold. However, if the emitting and 
Faraday-rotating regions are mixed, the Pi'R) distribution along 
the line of sight will no longer be a 5-function, and the re- 
constructed Pi'R) spectrum will be wider than the width of the 
RMSF. If the synchrotron-emitting and Faraday-rotating regions 
are fully mixed, and if the emissivity and the amount of Faraday 
rotation per parsec are independent of the position along the line 
of sight, then 

. (11) 



where "Rmax is the Faraday depth of the far side of the slab, 
and assuming that all radiation along the line of sight is emit- 
ted with a position angle of 0°. This result was first obtained by 
Burn ( 119661 ). and the distribution of synchrotron-emitting and 
Faraday-rotating regions along the line of sight that produces it 
has become known as a 'Burn slab'. Note that O again depends 
Unearly on A^, but that RM = 0.5'R^.^^. Sokoloff^ et al. ( fT998t 
have shown that a linear dependence of O on A^ occurs not 
only for a Burn slab, but for every distribution of synchrotron- 
emitting and Faraday-rotating regions that is symmetric along 
the line of sight between the observer and !R max. The recon- 
structed Pi'R) spectrum of a Burn slab has a finite extent, and 
may even be wider than the width of a single RMSF. 

It is important to realise that for a general distribution of 
emitting and Faraday-rotating regions, only those structures 
show up in the Pi'R) spectrum that illuminate a column of 
Faraday-rotating ISM from the back. If Faraday rotation occurs 
also at larger distances, but is not illuminated from the back by 
synchrotron emission, this will not appear in the PiK) spectrum. 

3. The data 

The present data set was obtained with the WSRT, a 14-element 
E-W interferometer of which 4 elements are moveable to im- 
prove (u,v) coverage. Each of the telescope dishes is 25 m in 
diameter The central coordinates of the region we study (in the 
constellation Gemini) are a = 7'' 18"' and 6 = 36°24' (J2000.0), 
which is Z x 181° and/? » 20° in Galactic coordinates. The 
GEMINI region was observed in 6 12 hour observing runs in 
December 2002 and January 2003 (see Table [1]). This yielded 
visibilities at baselines from 36 to 2760 meters, with an incre- 
ment of 12 meters. We tapered the individual frequency chan- 
nel maps in such a way that the synthesized beamsize for all 
maps is that of the 385 MHz beam of 2.76'x 4.70' (RA x DEC). 
Combining the 6 12 hr observing runs puts the first grating ring 
at 4. 1° (at 350 MHz) from the pointing centre, outside the 3°x3° 
area that we mapped for each individual pointing. 

We mapped an area of about 9°x9° with a 7x7 pointing mo- 
saic. In each night the same field was observed about 16 times, 
resulting in visibilities on 16 'spokes' in the (u,v) plane, each 
time integrating for 40 seconds before moving to the next field. 
The 1 .2° distance between pointing centres suppresses off-axis 
instrumental polarization to less than 1% (Wieringa et al. 1993] l. 
In our analysis we leave out the edges of the mosaic where in- 
strumental polarization effects are not suppressed by mosaick- 
ing. 

The observations cover the frequency range between 324 and 
387 MHz. 202 independent spectral channels are usable from 
a maximum of 224, each channel being 0.4 MHz wide (where 
we used a Hamming taper). 10 channels were flagged because 
Stokes V was contaminated by radio frequency interference, and 
one channel was manually flagged in all 6 12 hr runs. 

The data were reduced using the NEWSTAR data reduction 
package. Dipole gains and phases and leakage corrections were 
determined using the unpolarized calibrators 3C48, 3C147 and 
3C295. The flux scales of both unpolarized and polarized cali- 
brators are set by the calibrated flux of 3C286 (26.93 Jy at 325 
MHz - Baars et al. 11977 ). Due to an a-priori unknown phase off- 
set between the horizontal and vertical dipoles, signal can leak 
from Stokes U into Stokes V. We corrected for this by rotat- 
ing the polarization vector in the Stokes iU,V) plane back to the 
U axis, assuming that there is no signal in V. The polarized cali- 
brator sources 3C345 and DA240 defined the sense of derotation 
(i.e. to the positive or negative U-axis). Special care was taken to 
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Central position 
Mosaic size 
Pointings 
Frequencies 

Resolution 
Stokes V noise level 
K noise level^ 
Conversion Jy-K* 



(a. Show = (109° ,36.5"); (/, b) = (181° ,20") 

9°x 9° 
7x7 

324-387 MHz 

202 independent frequency channels (Hamming taper) 
2.76' X 4.70' 
6.2 mJy (2.0 K) 
0.5 mJy (0.14 K) 

1 mJy/beam = 0.32 K (at 345 MHz) 



Shortest baseline (m) 36 48 60 72 84 96 

Observing date 03/01/17 03/01/21 02/12/18 02/12/05 03/01/22 03/01/13 
(3T'/nim/dd) 

Start time (UT) 17:04 16:37 18:45 19:36 16:39 17:32 

End time (UT) 05:04 04:36 06:43 07:34 04:37 05:31 

Table 1. Characteristics of the GEMINI data set. Observing dates and times are given for each of the 12 hr observing runs, which 
have been indicated by their shortest baseline length. We calculated the conversion factor between mJy/beam and K at 345 MHz, 
the average of the sampling of the (usable) frequency channels in our data set. 




11+ 112 110 108 10^ 104 
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Fig. 3. Slices through the PCR) cube that show the strongest diffuse emission. The K of each sUce is indicated in rad/m^ in the upper 
right comer of each image. The RMSF has a FWHM of 12 rad/m^, and we sampled each RMSF with 2 images. The linear intensity 
scale is the same for all panels, and it saturates at 6.4 K. 



avoid automatic flagging of real signal on the shortest baselines. 
From Stokes V, which we assume to be empty, we estimate that 
the average noise level in the mosaics of the individual channels 
is 6.2 mJy (2.0 K). Instrumental polarization increases towards 



the edges of the maps, therefore we excluded these regions in 
determining the noise level. 

These observations were carried out in the evening and at 
night to Umit solar interference and to reduce the importance of 
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ionospheric RM variations. From 2 lines of sight with a strong 
polarized signal we estimate that the amounts of ionospheric 
Faraday rotation in the 6 nights are identical to within ~ 10°, 
so we did not correct for this. Ionosphere models indicate that 
the ionospheric contribution to RM during the observing nights 
was only 0.6 rad/m^ (Johnston-Hollitt, private communication). 

An interferometer will not cover all baseline lengths down 
to m, which means that maps of the sky that were made using 
an interferometer will miss structure on large angular scales. We 
will return to this point in Sect. 14.21 

4. Analysis 

4.1. The K datacube 

We calculated PCR) maps of the sky for Faraday depths from 
-1000 rad/m^ to +998 rad/m^ in steps of 6 rad/m-. As the 
width of the RMSF along the K axis of the datacube is about 
12 radians/m^ for our data set, this gives Nyquist sampling in 
Faraday depth. Due to the finite channelwidth of the data, the 
sensitivity of the PCR) spectra will drop at large K, as described 
by Eqn. [11 This becomes important for |'??| > 1.9 /SA^ ^ 1250 
rad/m^. 

In Fig. [3] we show slices through our datacube that show 
strong Galactic emission. All of the images saturate at 6.4 K. In 
Figs, m and |5] we summarise the information in the K datacube 
by plotting for each line of sight the maximum PCR) along that 
line of sight, and the !R at which this maximum occurs. We only 
plotted lines of sight where the main peak is more than twice as 
high as the second highest peak. From these figures it is clear that 
the majority of the lines of sight satisfies this condition. Because 
we treated the data in a slightly different way from what we de- 
scribed in Schnitzeler et al. ( |2007bl l, there are some minor dif- 
ferences with the corresponding figures in that paper From the 
PCR) maps at |'??| > 200 rad/m^ we estimate that the noise level 
in the PCR) slices is 0.5 mJy (0.14 K). 

The complexity of a line of sight depends on whether the 
main peak is resolved, and on how strong second and higher- 
order peaks are. The A criterion that we introduced in Schnitzeler 
et al. ( |2007bl l can be used to address the issue of resolution. We 
defined A as the root-mean-square vertical separation in the K 
spectrum between the main peak in the PCR) spectrum and the 
best-fitting RMSF, and we compare the observed A to the distri- 
bution of A that we simulated for an input signal of the appropri- 
ate strength + noise. In Fig.|6]we plot the lines of sight that have 
an unresolved main peak as green pixels. 'Unresolved' means 
that the A value of the main peak is < the A that we found for 
99% of our simulations of signal + noise the 3cr level of the 
Rayleigh distribution). 

Lines of sight with linear ^(A^) relations have special ge- 
ometries of Faraday-rotating and synchrotron-emitting regions, 
as discussed in Sect. |2l We therefore fitted a straight line to our 
<i>{A^) datapoints, where we included the periodicity of the data 
by using the procedure described in Schnitzeler et al. (2007b). 
Lines of sight with;^^^^^ < 2 are indicated in Fig.|6]as yellow pix- 
els, and lines of sight that have both a not measurably resolved 
main peak and a low;^^^^^ are indicated as red pixels. Since many 
lines of sight have a low;^'^^^^, but a high A, there are many yellow 
pixels, but relatively few red pixels. 

We simulated the PCR) response from the main peak plus 
noise by calculating the Fourier transform from Eqn. [3] for an 
input signal of a given strength + the noise levels that we derived 
from the Stokes V maps of the individual frequency channels. 
We assume that the Stokes V maps contain no signal. We ran our 
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Fig. 4. PCR) of the highest peak in the H spectrum. Each pixel 
represents one telescope beam, and the pixel size is 2.8'x4.9'. 
The lines of sight that are shown in this plot have a second peak 
in the PCR) spectrum that is at most half the strength of the main 
peak. Out of a total of 22.800 lines of sight, 2720 do not obey 
this selection criterion, and they are shown in black, as are the 
769 lines of sight that have PCR) < 2 K. The Icr noise level in 
the PCR) maps is 0.14 K. 
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Fig. 5. corresponding to the PCR) plotted in Fig.|4] The lines 
of sight plotted in these figures are selected by using the same 
criterion. Lines of sight indicated in white have K > 30 rad/m^. 
Lines of sight indicated in black either did not pass the selection 
criterion or have K < -6 rad/m^ (2720 and 727 lines of sight 
respectively). 



simulations for a grid of input signal strengths, that was matched 
to the strengths of the main peak in the PCR) spectra that we 
encountered in the data. In this way we calculated for each K in 
the spectrum the level below which 99% of the RMSF plus noise 
simulations lie. 
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Fig. 6. Distribution of independent lines of sight with a main 
peak that satisfies the A criterion (green pixels), or with;^^^^^ < 2 
and max(P('??)) > 2 K (yellow pixels), or that satisfy all three 
criteria (red pixels). Lines of sight that did not pass these criteria 
are shown as black pixels. The outer edges of the mosaic show 
an increase inx^^^, which explains the absence of yellow pixels. 



The main peak in the PCR) spectrum is not measurably re- 
solved for about 8% of the lines of sight in the current data set, 
in the sense that its value of A is lower than the A that we find for 
99% of our simulations of a RMSF + noise. This value is lower 
than the 14% quoted in Schnitzeler et al. ( I2007bl ). This differ- 
ence might be produced because we now use frequency channel 
maps that all have been tapered to the same beamsize. A not 
measurably resolved peak means that a Faraday-rotating region 
is physically separated from a region with synchrotron emission; 
mixing of the 2 regions produces broad peaks in the PCR) spec- 
trum. Since we only tested the main peak in our PCR) spectra, 
and no other features in the PCR) spectra, emission regions with 
a Faraday screen in front of them will be even more numerous 
than in 8% of the lines of sight. RFI and calibration errors can 
only increase A, which is another reason why there are probably 
more Faraday screens present in the ISM. 

Intuitively one would not expect such a separation between 
the emitting and Faraday-rotating regions along the line of sight: 
cosmic rays that produce synchrotron radiation are expected to 
be present wherever there is a magnetic field. Therefore one pos- 
sibility could be that we simply do not see the synchrotron emis- 
sion because the full magnetic field vector is pointing along the 
line of sight. The magnetic field component perpendicular to the 
line of sight, which determines the synchrotron emissivity, is in 
that case 0. Such a region would from our point of view then 
only produce Faraday rotation. Even though the strengths of the 
magnetic field component along the line of sight that we deter- 
mine in Sect. |5] are low, these are averages over the line of sight, 
and they are biased towards regions where the thermal (i.e. 10** 
K) electron density is higher 

The majority of lines of sight in our data have a secondary 
peak that is higher than what would be expected for one in lO'* 
noise realisations. However, the main peak is typically more than 
twice as high as the secondary peak, and thus dominates most of 
the lines of sight in the GEMINI region. This makes these lines 



of sight relatively simple, and we did not deconvolve the PCR) 
spectra. The reason why these lines of sight are not very complex 
could be because we are looking in the direction of the Galactic 
anti-centre, at an intermediate Galactic latitude. 



4.2. Missing large-scale structure 

Structure on large angular scales is missing from interferomet- 
ric observations if information on short baselines is not avail- 
able. Differences in Faraday rotation in between the source of the 
polarized emission and the observer can for neighbouring lines 
of sight 'transfer' polarized signal from large angular scales to 
small enough angular scales so that it can be picked up by the in- 
terferometer. Since this effect does not operate on Stokes /, this 
is the canonical explanation why structure in polarized intensity 
does not appear to have a counterpart in total intensity. In this 
section we expand on this idea, and we show that modulation by 
a linear gradient in Faraday depth "R shifts the angular frequency 
spectrum as a whole. 

For a linear variation of K with position / (i.e. a fixed gradi- 
ent VR), and wavelength A, this modulation of the polarization 
vector P can be expressed as 



Pobs(l) = PUl)e 



2iVRi-/ 



(12) 



where subscripts 'em' and 'obs' refer to the emission before and 
after it gets modulated by the foreground gradient in K. This lin- 
ear modulation with position shifts the angular frequency spec- 
trum as a whole: 



dl P,Ul) e-^""'' 

oo 

dl Pem(0 e2'^«^'''e-2-"' 

oo 



(13) 



f 



dl Pemil) e 



= FT{P,,J(u--V'RA^) 

7T 



(14) 



. One consequence of this is that foreground Faraday modula- 
tion makes the 0-angular frequency from /'em(0 visible if \VR\ 
is large enough. To make it detectable at the 25 m (projected) 
baseline length, which can be obtained with the WSRT, requires 
a |V!R| K 6 rad/m^/field of view at 350 MHz (where 1 field of 
view is 3° wide for our observations). Since we observe gradi- 
ents of this magnitude in the individual pointings of the GEMINI 
data set (Fig.|5j, Fig.|4]and Fig. |5] contain at least the 0-angular 
frequency = total polarized intensity component of /"em- 

The direction of the shift in Eqn. [14] is set by the sign of 
VR. Since the WSRT only measures one half of the (u,v)-plane, 
this means that the angular frequency spectrum could move 
away from the measurement points in the (u,v)-plane. However, 
^'obs(0 = 2obs(0 + if/obs(0, and since 2obs(0 and f/obs(0 are real 
quantities, their Fourier transforms are hermitian. The gradient 
from Eqn. [14] always moves both the positive and the negative 
angular frequencies of 2eni(0 and t/em(0 towards or away from 
the origin in the (u,v)-plane. Therefore, Pobsil) can be recon- 
structed even if only one half of the (u,v)-plane is observed. 

Note that foreground Faraday modulation shifts some an- 
gular frequencies from /"emO towards smaller (in an absolute 
sense) angular frequencies. At some point these angular frequen- 
cies from Pemil) will no longer be visible for an interferometer. 
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because it can only measure visibilities from a certain shortest 
baseline. For the analysis that we present in the remainder of 
this article it is sufficient that our data contain the total polar- 
ized intensity component of Pem- Single-dish observations are 
therefore still required to produce maps that contain all angular 
frequencies. 

If Faraday rotation in the foreground produces strong vari- 
ations in !R over the telescope beam, then this will also lead to 
beam depolarization. For a Gaussian beam, a gradient of more 
than 12 radians/m^/field of view is required to give more than 
10% beam depolarization (Sokoloff et al. 1 19981) , and such steep 
gradients are not observed in the GEMINI data. 

As a final note in this section, we showed in Schnitzeler et al. 
( |2007al l that the polarization angle gradients in the WENSS data 
in the region 130° </ < 170°, -5° < b < 30° ai-e on average about 
2 radians/m^/° along Galactic latitude, which translates into the 6 
radians/m^/field of view required to shift the angular frequency 
spectrum towards observable frequencies. These WENSS gra- 
dients therefore also show at least the total polarized intensity 
component of the angular frequency spectrum. 



5. The line-of-sight component of the large-scale 
magnetic field 

In this section we derive the strength of the magnetic field com- 
ponent along the line of sight, B^. We use the Faraday depths 
from Fig. |5] and we determine the contribution from thermal 
electrons to (Eqn. |2]l by converting Ha intensities from the 
WHAM survey (Haffner et al. 120031) first to emission mea- 
sures EM - nidi, and then to dispersion measures DM - 

J "CO 
^ riedl, using an empirical relation between EM and DM that 

Berkhuijsen et al. ('2006') derived for Galactic pulsars. The inte- 
grals for EM and DM are over the entire line of sight, [EM] = 
cm"''pc, [DM] - cm^-'pc, [rig] = cm"-' and [d/] - pc. 

To estimate DM, we could also have used the electron den- 
sity models by Reynolds ( 119911 ) or by Cordes & Lazio ( 120031 ) as 
alternatives to the relation between EM and DM that Berkhuijsen 
et al. derived. In the Reynolds model, 40% of the line of sight 
contains free electrons, at a density of 0.08 cm"^, and there are 
no free electrons in the remainder of the line of sight. Cordes 
& Lazio fitted a smooth electron density model consisting of a 
thin and thick disc + spiral arms to the DM measured for about 
1200 Galactic pulsars and 100 extragalactic sources. Small-scale 
structure in the Galactic ISM is included in their model as a list 
of local clumps and voids. Berkhuijsen et al. used DM measured 
for a sample of 157 pulsars, and accurately determined the EM 
between us and the pulsar by correcting for the EM that is pro- 
duced beyond the pulsar, and for reddening occurring in front 
of the pulsar By combining Berkhuijsen et al.'s statistical de- 
scription of the ISM with the WHAM Ha intensities that are 
measured on as; l°xl° grid, we can more accurately determine 
DM in specific directions than is possible with the models by 
Reynolds or by Cordes & Lazio. 

Significant numbers of free electrons are present in both the 
warm ionized medium (WIM) and hot intercloud medium (HIM) 
phases of the ISM. Snowden et al. ( 1997 ) modelled the Galactic 
X-ray halo from ROSAT diffuse X-ray background maps. Their 
model consists of a cylindrical bulge of radius 5.6 kpc, with a 
HIM electron density that decays exponentially away from the 
Galactic midplane. The HIM electron density in their model is 
in the Galactic midplane a factor of 10 smaller than the aver- 
age WIM electron density (3.5 x 10""^ cm"-* compared to ~ 3 
xlO"^ cm"-*). If we assume that the HIM electron densities in 



the outer Galaxy (in the direction of GEMINI) are comparable 
to the densities in the inner Galaxy, then the DM produced in the 
HIM is small compared to the DM of the WIM, and we therefore 
neglected the former 

In the top panel of Fig. [T] we plot the Ha intensities (in 
Rayleigh) for WHAM lines of sight in the GEMINI region. To 
convert the Ha intensities to EM, we use Eqn. 1 from Haffner et 
al. (1998): 



EM = 2.757^ '' iHa e 



2.2 E(B-V) 



(15) 



. T4 is the temperature of the WIM gas in units of 10"* K, 
which is typically 0.8 (Reynolds 1985). Ihu is the Ha intensity 
in Rayleigh, and E(B-V) is the interstellar B-V reddening. We 
determined EM for each of the WHAM Unes of sight by using 
reddening values from Schlegel et al. d 19981 1 and the programs 
available on Marc Davis' websit^H. 

The 157 pulsars with \b\ > 5° in the sample by Berkhuijsen 
et al. follow the relation 



EM = (0.042 + 0.014) DM 



1.47±0.09 



(16) 



. With the K from Fig. |5] and the DM that we calculated from 
the WHAM Ha intensities, we can then calculate the electron- 
density weighted average By along the line of sight, (By): 



<B||> ^ 



/ rieB ■ dl 

line of sight K 



j riAl 0.81 DM 

line of sight 



(17) 



where is in fiG. In the bottom panel of Fig. |7] we show 
the that we calculated in this way. Since the WHAM beam 
(FWHM a! 1°) is much larger than the synthesized WSRT beam, 
we used the P('??)-weighted H average over the WHAM beam to 
calculate (B||). Also note that since Schlegel et al. give redden- 
ing values for the entire line of sight through the Galaxy, the EM 
and DM that we calculate are also values for the entire line of 
sight. Since we do not know where the polarized emission that 
we observe originates, the KB||)| that we derive from Eqn. [Tt] are 
lower limits to the actual KB||)|. 

An upper limit for B\\ can be estimated by using a pitch an- 
gle of 8° and a strength of 4 //G for the large-scale magnetic 
field (both values from Beck '2007a'), which gives Breg,||,max ~ 0.6 
jjG. In the presence of a random magnetic field, By approaches 
the Bieg.ti component of the regular field if either the reversals in 
the direction of the random field are unimportant, or if the line 
of sight is long enough so that the contributions to H from such 
small-scale reversals cancel each other along the line of sight. As 
discussed in Haverkorn et al. (|2004| l, the random magnetic field 
has a correlation length of order 10 parsec. Shortly we will esti- 
mate that the length of the line of sight through ionized regions 
in the direction of GEMINI is of the order of 130 parsec. The 
number of 'draws' of the random field is therefore small, which 
means that the random magnetic field will have a non-negligible 
contribution (compared to the large-scale field) to the {B\\) that 
we derive. The (By) in the bottom panel of Fig. |7] already show 
that the (B|j) in GEMINI vary much more than what we would 
expect on the basis of the large-scale magnetic field alone. 

The (By) in the bottom panel of Fig. |7]are smaller than or 
comparable to the Breg,y,max that we estimated. The (By) we de- 
rive are small, which can be expected since we are looking in the 
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Fig. 7. WHAM Ho- intensities in Rayleigh (top panel) and the 
electron-density weighted strengths of the magnetic field com- 
ponent parallel to the line of sight (By) [fiG] that we derived 
(bottom panel). Open circles indicate positive there are no 
negative (By) < -0.1 juG in this figure. KB||)| < 0.1 fiG are shown 
as '+'. 

direction of the Galactic anti-centre, where the regular magnetic 
field is mostly perpendicular to the line of sight. It is clear from 
the bottom panel in Fig. |2]that variations in EM alone are not 
enough to explain the measured variation in but that varia- 
tions in Z?|| are required as well. Further support for this is given 
by the fact that K in Fig. |5] changes sign in the area covered by 
the GEMINI data, which can only be explained by a change in 
the direction of By. The By structure in Fig.|5]is not associated 
with any of the 3 largest non-thermal loops in the Galaxy, the 
North Polar Spur, Cetus Arc and Loop III. The locations of these 
structures are described in Elliott ( 1970]l. 

Berkhuijsen et al. also constructed a model of the electron 
density distribution along the line of sight, that we can use to es- 
timate the length of the line of sight that produces the observed 
EM. Similar to the model by Reynolds ( 119911 ), their model con- 
sists of cells that are empty or filled with a constant electron den- 
sity. Contrary to the Reynolds model, the electron density of the 
filled cells ('clouds'), We.cioud, is allowed to vary between lines of 
sight: 



«e,cloud= ^ =0.12EM°-^2 (j8) 

DM 

where we used Eqn. [16] Berkhuijsen et al. found that the fill- 
ing factor / of electron clouds depends inversely on the electron 
density in the clouds: / = (0.0184 + 0.001 1) n;l^f-°\ We can 
use this relation to determine the length of the line of sight D (in 
parsec) that produces the observed EM; note that D includes both 
the fraction of the line of sight that passes through the electron 
clouds and also the electron-free regions between the clouds: 



D 



DM^ 



/(em 




404 EM 



0.70 



(19) 



10 108 
RA(2000) ["] 

Fig. 8. RM distribution for the polarized point sources from ta- 
ble|2] The size of the circles is proportional to RM, shown in the 
scale on top of the figure in units of radians/m^. Black circles 
indicate negative RM, and white circles indicate positive RM. 
|RM| < 4 rad/m^ are missing from this figure because of our 
selection criteria. We excluded the grey region when looking for 
polarized point sources, because instrumental polarization levels 
are higher there. 



. E(B-V) = 0.071 magnitudes for (a, 6) = (110°,37°), which 
means that the average Ha intensity in GEMINI of 1.78 R trans- 
lates into an EM = 4.68 cm^^pc. This gives ne.cioud - 0.19 cm"-', 
/ = 0. 1 1 , and D = 1 .2 kpc. This filling factor and D imply that the 
combined pathlength through the electron clouds is fxD x 130 
parsec. We calculated EM and DM for the entire line of sight, 
but it is unclear whether we can observe the polarized emission 
coming from the furthest parts along the line of sight. D is there- 
fore an upper limit to the length of the line of sight that produces 
the!RinFig.|5] 



6. Polarized point sources 

Since the line of sight towards extragalactic sources passes 
through the entire Milky Way, whereas diff'erent depolarization 
effects can limit the line of sight of the diffuse emission, the 
information we obtain from the of extragalactic sources can 
complement what we learn from the diffuse emission. However, 
extragalactic sources will also have an intrinsic fi. 

To find polarized point sources, we first made maps that only 
contain baselines > 250 m, to filter out the strong diffuse emis- 
sion from our PCR) datacube. These maps use a Gaussian taper 
that reaches a value of 0.25 at a baseline length of 2500 m. We 
then Fourier transformed these channel maps from space to K 
space, and constructed PCR) maps from = -1000 rad/m^ to H 
- +998 rad/m^, in a way identical to what we described in Sect. 

KB 

The polarized intensity of an extragalactic source with rota- 
tion measure RM will be attenuated by bandwidth depolarization 
according to sinc(|RM|5/l^), if the polarization angle varies lin- 
early with A^. Here 6A^ is the A^ width of the frequency channels. 
The channels that we use have on average a 6A^ ^ 1.5 x lO"-' 
m^, therefore bandwidth depolarization of a point source with 
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Table 2. Properties of the polarized point sources we found in 
our data. Shown are the equatorial coordinates (J2000.0) of the 
source in decimal degrees, the RM that we fitted and its error 
(both in rad/m^), the reduced ^xl^A^ of the fit, the maximum 
PiU) along the line of sight and the total intensity / (both in 
mJy), and the polarization fraction (in %). 



\'R\ - 1000 rad/m^ will reduce the signal to 66%, which would 
still be detectable. 

We looked for each line of sight in the PCR) datacube if 
1) there was a Stokes / counterpart brighter than 150 mJy, 2) 
if the maximum Pili) along the line of sight was more than a 
5cr detection and 3) if the ratio of the maximum f (!R) along the 
line of sight/ total intensity / was larger than 0.02, meaning that 
this source is more than 2% polarized (the instrumental polar- 
ization level is about 1%). We also exclude sources with \'R\ < 4 
rad/m^ because instrumentally polarized point sources will show 
up in the PCR) datacube at rad/m^, and the finite width of the 
RMSF along the PCR) axis will produce a strong signal also in 
the vicinity of rad/m^. Furthermore we excluded a region of 
about 0.9° from the edge of the mosaic, where there is no over- 
lap between pointings, and instrumental polarization levels are 
therefore much higher there. We then fitted a RM to the <1>(/1^) 
distribution of the sources that satisfy these criteria in the way 
described in Schnitzeler et al. ( |2007b| . 

In Table|2]we list the properties of the sources that we found, 
and we plot the RM of these sources in Fig. [8] The low x^^^ of 
the RM fits indicate that there is only one PCR) component along 
the K axis, and that we have effectively filtered out the diffuse 
Galactic contribution to the PCR) spectrum. This, in combina- 
tion with the fact that these sources are detected in maps where 
we left out low angular frequencies on the sky, means that the 
sources that we detected are compact in all 3 dimensions of the 
PCR) cube, and are therefore likely to be either polarized pulsars 
or polarized extragalactic sources. Since the surface density of 
the latter is so much higher than that of the former, the sources 
that we detected are probably extragalactic in origin. 

Most of the sources in the upper left quadrant of Fig. |8] show 
positive RMs of x 5-15 rad/m^, which is similar to the Faraday 



depths of the diffuse emission from Fig. |5] In other regions the 
source density is insufficient to reach definite conclusions. The 
RM that we observe for an extragalactic source is the combina- 
tion of a Galactic RM and a RM that is intrinsic to the source. 
However, since a number of sources show the same RM, and 
since the intrinsic RM of these sources are uncorrected, what 
we see is probably mostly the Galactic RM contribution for the 
entire line of sight through the Galaxy. 



7. Depolarization canals in PCR) maps 

Dark, narrow channels ('canals') can be clearly seen in the im- 
ages of Fig. |3] Canals have been known to exist in P(A) images 
(e.g. Haverkorn et al. 2000), but, as far as we know, canals in 
PCR) maps have only been mentioned by De Bruyn et al. (2006]l. 

In Fig.|9]we show slices through the PCR) cube that have par- 
ticularly strong depolarization canals. It appears from this figure 
that the canals essentially are not changing position with chang- 
ing K. This is demonstrated by the reference crosses that we po- 
sitioned on three canals in this figure. In the first and second rows 
of Fig. [TO]we plot the P(A^) (left column) and <S>(A^) (right col- 
umn) spectra for the two canal pixels from the top row of Fig|9] 
('CH r and 'CH 2' respectively). In the bottom row of Fig. [10] 
we plot the spectra for the canal pixel in the bottom row of Fig. 

m 

In Fig.[TO]we also plot the P{A^) and (i>{A^) spectra for 2 lines 
of sight that lie on either side of the central canal pixel, at half 
a (synthesized) beamwidth away. These adjacent lines of sight 
share pixels with the line of sight that is centred on the canal 
pixel, and with these lines of sight we can investigate the nature 
of the canal pixels. 

The O(A^) spectra of the neighbouring lines of sight were 
unwrapped by fitting a straight line to the polarization angles 
by hand. As a first guess for the slope of this line, we used 
the "R at which the PCR) of the neighbouring lines of sight is 
high. Furthermore we required that for each line of sight the dif- 
ference in polarization angle between the successive (plotted) 
wavelength^ channels must be smaller than 90°. 

The difference in (^(A^) for lines of sight adjacent to the 
canal pixels is large, often even almost 90°, averaged over the 
available A^. Also, the polarized brightness temperature of the 
2 neighbouring lines of sight of each canal pixel are on average 
comparable. These 2 facts combined indicates that beam depo- 
larization, i.e. summing polarization vectors with similar lengths 
but different orientations over the (synthesized) telescope beam, 
is an important factor in creating canals in the PCR) maps. In 
Haverkorn et al. ( 2000) the same mechanism was identified to 
be producing canals in PiA^) maps. Line-of-sight depolarization, 
occurring for example as the nulls in Eqn. [TT] produces canals 
only at certain wavelength^. Shukurov & Berkhuijsen (2003]l 
suggested that this mechanism is producing canals that are seen 
towards M31. This type of canal is characterised as extended 
along the "R axis of the PCR) cube, contrary to the properties of 
the canal pixels we present here. 



8. The area around {l,b) = (109°,34.5°) 

The 2°x2° area centred on (l,b) = (109.5°, 34.5°) is bright in 
polarized intensity (Fig.|4]i, shows a very uniform distribution in 
K (Fig. |5j, and has a relatively small surface density of strong 
depolarization canals. Furthermore the x^^^ fitted to the <i>{A^) 
spectra are low (Fig. |6]i: 30% of the lines of sight in this area 
have a , < 2, and nearly 75% have x^- a ^ ^- There is no 



D.H.F.M. Schnitzeler et al: WSRT Faraday tomography of the Galactic ISM at A ~ 0.86 m 



11 



40.0 - 



39.5 - 



o 

o 

o 

39.0- 

LjJ 

Q 



38.5 - 



^4 




1 1 2.0 1 1 1 .5 1 1 1 .0 
RA(2000) [°] 




1 1 2.0 1 1 1 .5 1 1 1 .0 
RA(2000) [°] 

Fig. 9. Depolarization canals in the PCR) datacube. The K [rad/m^] of each of the PCR) slices is indicated in the lower left. All maps 
have the same linear greyscale that saturates at 4.8 K. The area shown in the top row of this figure lies in the upper left of Fig. [3] 
For the 2 pixels that are indicated by crosses, and that lie in depolarization canals, we plot the P{A^) and <i>(A^) spectra in the first 
row of Fig. [To] (left cross, 'CH 1' in Fig. [TOt and middle row of Fig. [10] (right cross, 'CH 2'). The bottom row shows a canal in a 
different part of the mosaic, and the P{A^) and <^(A^) spectra of the canal pixel that is indicated by a cross are shown in the bottom 
row ofFig.[T0]('CH3'). 



enhanced Ha emission in this region (see the top panel of Fig. 
[7]), and B|| varies only little (bottom panel of Fig.|7]i. The single- 
dish 408 MHz total intensity map by Haslam et al. ( 1982) shows 
hardly an increase in brightness temperature in this region. There 
are not enough lines of sight in the Brouw & Spoelstra ( 1 19761 1 
data set to discern the area around il,b) = (109.5°, 34.5°) from 
its surroundings. 

The increased polarized brightness temperature with respect 
to its surroundings could be due to enhanced synchrotron emis- 
sion, or a local decrease in depolarization, or to a combination 
of these. The polarized brightness temperature of this region is 
more than 3 K larger than that of its surroundings, and this would 
translate into a brightness temperature difference of at least 3 
K/0.7 - 4.3 K, assuming that the radiation is 70% polarized. 
Lower polarization percentages increase this brightness tem- 
perature contrast even further An enhancement in synchrotron 
emission can be ruled out as an explanation for the increased 
T^b.poi of this area, since this is not seen in the Haslam data, and 
since beam averaging is not strong enough to attenuate the total 
intensity signal (the Haslam beam measures 0.85°x0.85°). 

The small variation in K over this region, combined with 
an underabundance of strong depolarization canals, furthermore 



indicates that depolarization across the synthesized telescope 
beam is not a strong effect. This leaves a decrease in the amount 
of depolarization along the line of sight as an explanation for the 
observed increase in Tb.poi. The low;^'^^^ of lines of sight in this 
area, in combination with secondary peaks that are not stronger 
than half the strength of the main peak, also indicates that the 
lines of sight in this area are not very complex, and line-of-sight 
depolarization effects are therefore not very important. We thus 
think that the origin of this particular region lies in a decreased 
amount of depolarization compared to its surroundings. 

9. Conclusions 

We applied Faraday tomography to high spectral-resolution ra- 
dio polarization data to study the properties of the magnetised 
Galactic ISM. We showed that differential Faraday modulation 
in the foreground can shift the emitted angular frequency spec- 
trum as a whole. In particular the 0-angular frequency from the 
emitted radiation is modulated in this way towards smaller angu- 
lar scales, that are observable with an interferometer We quan- 
tified how strong a linear gradient in fl should be to accomplish 
this, and we showed that the gradients in that we find for the 
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Fig. 10. Polarized brightness temperature Tb,poi [K] and polarization angle O [degrees] spectra for 3 canal pixels, and two adjacent 
pixels half a beamwidth from each canal pixel. Top row: canal pixel at (a, (5)=(111.24°,39.72°), and two adjacent lines of sight at 
□'=111.27° (dotted line), and q'=1 11.21° (dashed line). Middle row: canal pixel at (a, 5)=(1 10.98°,39.57°), and adjacent lines of 
sight at Qr=111.0r dotted line, and a'=110.95°, dashed line. Bottom row: canal pixel at (a;,5)=(111.50°,34.94°), and two adjacent 
lines of sight at q'=1 1 1.53°, dotted line, and a- 1 1 1.47°, dashed line. The Icr errors for polarized intensity and polarization angle 
are 2 K and 15° (for a S/N = 2 detection) respectively Every 8* channel is plotted. Bad frequency channels were selected using the 
same selection criteria as in Sect. |3] and discarded. 



main peak in the PCR) spectra are indeed strong enough so that 
we are sensitive to the 0-angular frequency that is emitted at 
these Faraday depths. 



The main peak in the PCR) spectrum is not measurably re- 
solved for 8% of the lines of sight in our data set. An unresolved 
peak means that there is no synchrotron emission occurring in 
the Faraday-rotating region. This is unexpected, since cosmic 
rays that produce synchrotron emission pervade the ISM. We 
propose that the magnetic field orientation plays an important 
role in this, in that a magnetic field that is oriented along the 
line of sight produces only synchrotron emission perpendicular 
to the line of sight, hiding it from our view. By using the ob- 
served emission measures from the WHAM survey as input, and 
by simulating the thermal electron contribution to ^ , we could 
determine the magnetic field strength along the line of sight, and 
we mapped this quantity over the area covered by our GEMINI 
data. The polarized point sources we found in our data have 
RMs that are comparable to the we find for the diff'use emis- 
sion. In our PCR) maps we found depolarization canals, narrow 
structures with only a small fraction of the polarized intensity 
of adjacent lines of sight. We established that for a number of 
deep canals the polarization angle changes in many frequency 
channels by a large amount over the canal, in some cases 90°. 
We therefore argue that depolarization over the synthesized tele- 
scope beam is producing (at least) these canals. Finally, we in- 
vestigated the properties of a large, conspicuous area in our data, 
and we argue that a decrease in depolarization along the line of 
sight as compared to its surroundings is probably responsible for 
the observational features of this region. 
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